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This article describes three Mathematics packages for the automatic calculation of one-loop Feynman diagrams: 
the diagrams are generated with FeynArts, algebraically simplified with FormCalc, and finally evaluated numer- 
ically using the LoopTools package. The calculations are performed analytically as far as possible, with results 
given in a form well suited for numerical evaluation. The latter is straightforward with the utility programs 
provided by FormCalc (e.g. for translation into Fortran code) and the implementations of the one-loop integrals 
in LoopTools. The programs are also equipped for calculations in supersymmetric models. 



1. Introduction 



The precision of experimental data achieved at 
present colliders has in many cases reached or ex- 
ceeded the per cent level. Obviously a compara- 
ble accuracy on the theoretical side is needed in 
order to draw significant conclusions from such 
precise measurements. For many observables this 
means that a one-loop calculation is the lowest 
acceptable approximation. 

Doing one-loop calculations by hand is labo- 
rious and error-prone and in some cases simply 
impossible. So for some time already, software 
packages have been developed to automate these 
calculations (e.g. [ |], ||). Incidentally, full au- 
tomation is possible only up to one loop since 
no algorithms generic enough for the computa- 
tion of arbitrary multi-loop Feynman diagrams 
are known at present. One remaining obstacle is 
that the existing packages generally tackle only 
part of the problem, and one still has to spend 
considerable effort adapting conventions etc. to 
make them work together. 

In this paper the three Mathematica pack- 
ages FeynArts, FormCalc, and LoopTools are pre- 
sented which work hand in hand. The user has 
to supply only small driver programs whose main 
purpose is to specify the necessary input parame- 
ters. This makes the whole system very "open" in 
the sense that the results are returned as Math- 
ematica expressions which can easily be manipu- 
lated, e.g. to select or modify terms. 
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Figure 1 . Steps in a one- loop calculation and the 
distribution of tasks among the programs Feyn- 
Arts, FormCalc, and LoopTools. 
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Since one-loop calculations can range anywhere 
from a handful to several hundreds of diagrams 
(particularly so in models with many particles 
like the MSSM), speed is an issue, too. Form- 
Calc, the program which does the algebraic sim- 
plification, therefore uses FORM for the time- 
consuming parts of the calculation. Owing to 
FORM's speed, FormCalc can process, for ex- 
ample, the 1000-odd one-loop diagrams of W-W 
scattering in the Standard Model [ |(| in about 5 
minutes on an ordinary Pentium PC. 

The following sections describe the main func- 
tions of each program. Furthermore, the Form- 
Calc package contains two sample calculations in 
the electroweak Standard Model, ZZ -> ZZ [ § 
and e + e~ — > tt [ [|, which demonstrate how the 
programs are used together. 

2. FeynArts 

FeynArts is a Mathematica package for the 
generation and visualization of Feynman dia- 
grams and amplitudes. The current version 2.2 is 
a much-expanded version of FeynArts 1 [ |9| . The 
two most important new features are the genera- 
tion of counter-term diagrams and the ability to 
deal with supersymmetric theories (cf. Sect. |5|). 

FeynArts works in three basic steps, sketched 
in Fig. g. The first step is to create all different 
topologies for a given number of loops and ex- 
ternal legs. For example, to create all one-loop 
topologies for a 1 — > 2 process, the following call 
to CreateTopologies is used: 

top = CreateTopologies [1 , 1 -> 2] 

In the second step, the actual particles in the 
model have to be distributed over the topologies 
in all allowed ways. E.g. the diagrams for Z — > bb 
are produced with 

ins = InsertFields [ top, 

V[2] -> {F[4,{3}], -F[4,{3}]> ] 

where F[4,{3}] is the 6-quark, -F[4,{3>] its 
antiparticle, and V [2] the Z boson. 

The fields, their propagators, and their cou- 
plings are defined in a special file, the model 
file, which the user can supply or modify. The 
following model files are included in FeynArts: 
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(CreateTopologies) 
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Figure 2. Flowchart for the generation of Feyn- 
man amplitudes with FeynArts. 



the electroweak Standard Model (SM.mod) [Q, 
the same including QCD (SMQCD . mod) , and in 
the background-field formulation (SMbgf . mod) . 
These model files all include the full set of one- 
loop counter terms. A model file for the MSSM 
is in preparation. 

The diagrams can be drawn with Paint [ins] , 
depending on the options either on screen, or in 
a PostScript or FTgX file. 

Finally, the analytic expressions for the dia- 
grams are obtained by 

amp = CreateFeynAmp [ins] 

An important feature of FeynArts is that it dis- 
tinguishes three levels of fields: 

- Generic level, e.g. the fermion F, 

- Classes level, e.g. the down-type quark F [4] , 

- Particles level, e.g. the 6-quark F[4, {3}] . 
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This is useful for two reasons: 

The kinematic structure of a coupling is fixed 
once the generic fields are specified. For exam- 
ple, all fermion-fermion-scalar couplings are of 
the form 



C(F, F, S) = G- 



where u>± = (1 ± 7s)/2 are the chirality projec- 
tors. This means that most algebraic simplifica- 
tions like the tensor reduction need to be carried 
out on the Generic-level amplitude only. 

Furthermore, it is more economic to perform in- 
dex summations (e.g. over the fermion-generation 
index) in a loop over Classes-level amplitudes 
instead of generating many Particles-level dia- 
grams. 

3. FormCalc 

FeynArts produces very symbolic output which 
cannot straightforwardly be implemented in a 
numerical program. Its evaluation proceeds in- 
stead in two steps: first, algebraic simplification 
in Mathematics; then, translation into a Fortran 
program which computes the squared matrix ele- 
ment and from this the desired quantities (cross- 
sections, decay rates, asymmetries, etc.). 

3.1. Algebraic simplification 

The symbolic expressions for the diagrams are 
simplified algebraically with FormCalc which re- 
turns the results in a form well suited for numer- 
ical evaluation. Specifically, FormCalc performs 
the following simplifications: 

- indices are contracted as far as possible, 

- fermion traces are evaluated, 

- open fermion chains are simplified using the 
Dirac equation, 

- colour structures are simplified using the 
SU(iV) algebra, 

- the tensor reduction is done, 

- the results are partially factored, 

- abbreviations are introduced. 

The internal structure of FormCalc is simple: it 
prepares the symbolic expressions of the diagrams 
in an input file for FORM, runs FORM, and re- 
trieves the results via the MathLink interface (see 
Fig. |J). This is done completely without inter- 
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Figure 3. The interplay between Mathematica 
and FORM in FormCalc. 



vention by the user, i.e. the user does not see 
the FORM code. FormCalc thus combines the 
speed of FORM with the powerful instruction set 
of Mathematica and the latter greatly facilitates 
further processing of the results. 

The main function in FormCalc is OneLoop (the 
name is not strictly correct since it works also 
with tree graphs). It is used like this: 

« FormCalc' 

amps = << myamps.m 

result = OneLoop [amps] 

where it is assumed that the file myamps . m con- 
tains amplitudes generated by FeynArts. Note 
that OneLoop needs no declarations of the kine- 
matics of the underlying process; it uses the in- 
formation FeynArts hands down. 

Even more comprehensive than OneLoop, the 
function ProcessFile can process entire files. It 
collects the diagrams into blocks such that in- 
dex summations (e.g. over fermion generations) 
can later be carried out easily, i.e. only diagrams 
which are summed over the same indices are put 
in one block. This nicely complements the genera- 
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tion of Classes-level diagrams in FeynArts, which 
leaves the index summations to the numerical 
evaluation in order to reduce the number of di- 
agrams. ProcessFile is invoked e.g. as 

ProcessFile ["vertex . amp" , "vertex"] 

which reads the FeynArts amplitudes from the 
input file vertex . amp and produces output files of 
the form vertexid.m, where id is some identifier 
for a particular block. 

The output of OneLoop or ProcessFile is in 
general a linear combination of loop integrals with 
prefactors that contain model parameters, kine- 
matic variables, and abbreviations introduced by 
FormCalc. Such abbreviations are introduced for 
spinor chains, scalar products of vectors, and ep- 
silon tensors contracted with vectors. A term in 
the output could for instance look like 

C0i[cc0, MW2, MW2, S, MW2, MZ2, MW2] * 
( -4 Alfa2 CW2 MW2/SW2 S AbbSuml6 + 
32 Alfa2 CW2/SW2 S 2 AbbSum28 + 
4 Alfa2 CW2/SW2 S 2 AbbSum30 - 
8 Alfa2 CW2/SW2 S 2 AbbSum7 + 
Alfa2 CW2/SW2 S(T-U) Abbl + 
8 Alfa2 CW2/SW2 S(T-U) AbbSum29 ) 

The first line represents the one-loop integral 
C (M^,M^,s,M^,M^,M^), which is multi- 
plied with a linear combination of abbreviations 
like Abbl or AbbSum29 with certain coefficients. 
These coefficients contain kincmatical variables 
like the Mandelstam variables S, T, and U and 
model parameters, e.g. Alf a2 = a 2 . 

The automatic introduction of abbreviations is 
a very important feature which can drastically re- 
duce the size of an amplitude, particularly so be- 
cause the abbreviations are nested in three levels. 
Here is an example: 



AbbSum29 = Abb2 + Abb22 + Abb23 + Abb3 



Abb22 = Pairl Pair3 Pair6 



|Pair3 = Pair [e [3] , k [1] ] 

The definitions of the abbreviations can be re- 
trieved by Abbreviations [] which returns a list 
of rules such that result //. Abbreviations [] 
gives the full, unabbreviated expression. 



3.2. Translation to Fortran code 

For numerical evaluation, the Mathematica ex- 
pressions produced by FormCalc need to be trans- 
lated into a Fortran program. (They could, in 
principle, be evaluated in Mathematica directly, 
but this becomes rather slow for large ampli- 
tudes.) The translation is done by the program 
NumPrep, which is part of the FormCalc pack- 
age. The philosophy of NumPrep is that the user 
should not have to modify the generated code. 
This means that the code has to be encapsulated 
(i.e. no loose ends the user has to bother with), 
and that all necessary subsidiary files (include 
files, makefile) have to be produced, too. 

From the point of view of the Fortran program- 
mer who wants to use the generated code in his 
program, the output of NumPrep is a single sub- 
routine called 



squared_me (fci, 



, fcjv, £!)•••, £n, Ai, 



which takes as input the external momenta, po- 
larization vectors, and helicities, and returns the 
squared matrix element. To obtain actual numer- 
ical results from the generated code, one needs in 
addition a driver program whose task is to initial- 
ize the model parameters, set up the kinematics, 
invoke the squared_me subroutine, perform nec- 
essary phase-space integrations, and finally write 
the results to a file. A sample driver program for 
2^2 processes is included in FormCalc. 

Finally, the generated code has to be linked 
with the LoopTools library which provides the 
one-loop functions. 

4. LoopTools 

LoopTools supplies the actual numerical imple- 
mentations of the one-loop integrals needed for 
programs made from the FormCalc output. It is 
based on the reliable package FF [ [jij and pro- 
vides in addition to the scalar integrals of FF also 
the tensor coefficients in the conventions of [ [ll[] . 
LoopTools offers three interfaces: Fortran, C++, 
and Mathematica, so most programming tastes 
should be served. 

Using the LoopTools functions in Fortran and 
CH — h is very similar. In Fortran it is necessary 
to include the file looptools.h in every func- 
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tion or subroutine (for the common blocks). In 
CH — hj clooptools .h must be included once. Be- 
fore using any LoopTools function, ff ini must 
be called and at the end of the calculation f f exi 
may be called to obtain a summary of errors. It is 
of course possible to change parameters like the 
scale fi from dimensional regularization; this is 
described in detail in the manual [ [l^] . 

A very simple Fortran program would for in- 
stance be 

program simple 
#include "looptools .h" 
call ffini 

print *, B0(1000D0,50D0,80D0) 

call ffexi 

end 

The C++- version of this program is 
#include "clooptools .h" 

mainO 
{ 

ffini () ; 

cout « B0(1000. ,50. ,80.) « "\n"; 
ffexiO ; 

> 

The Mathematica interface is even simpler to use: 
In[l]:= Install ["LoopTools"] 
In [2] := BO [1000, 50, 80] 
Out [2]= -4.40593 + 2.70414 I 

5. Calculations in Supersymmetric Models 

Special emphasis has been placed on the possi- 
bility to do calculations in supersymmetric mod- 
els with FeynArts and FormCalc. In particular 
the following two fundamental problems become 
relevant in supersymmetric theories: 

Problem 1: SUSY theories in general contain 
Majorana fermions and hence fermion-number- 
violating couplings (e.g. quark-squark-gluino) . 
The textbook prescription of ordering the Dirac 
matrices opposite to their occurrence along the 
arrows on fermionic lines obviously breaks down 



in this case since one cannot define a fermion- 
number flow. (Put differently, Majorana-fcrmion 
lines have no arrow.) 

Solution: FeynArts uses the "flipping-rule" 
algorithm [ |l3| : instead of traversing the fermion 
lines along the fermion-number flow imposed from 
the outside, it chooses a particular direction for 
each fermion chain. If it later turns out that, for 
a Dirac fermion, the chosen direction is opposite 
to the actual fermion flow, it applies a so-called 
flipping rule. 

Problem 2: Dimensional regularization, the 
default regularization scheme employed by Form- 
Calc, is known to break supersymmetry [ [l4| . 

Solution: FormCalc has two regularization 
schemes built in which are chosen with the vari- 
able $Dimension. The default is $Dimension = 
D which corresponds to dimensional regulariza- 
tion. Putting $Dimension = 4 switches to con- 
strained differential renormalization [ ^| . The lat- 
ter technique is equivalent at the one-loop level to 
regularization by dimensional reduction [ Q and 
is hence suited for calculations in SUSY models. 

6. Requirements and Availability 

All three packages require Mathematica 2.2 
or above; FormCalc needs in addition FORM, 
preferably version 2 or above; LoopTools needs 
a Fortran compiler and gcc/g++. 

The packages should compile and run without 
change on any Unix platform. They are specifi- 
cally known to work under DEC Unix, HP-UX, 
Linux, Solaris, and AIX. A comprehensive man- 
ual which explains installation and usage is in- 
cluded in each package. All three packages are 
open-source programs and stand under the GNU 
library general public license. They are available 
from 



http 


//www 


f eynarts 


de 




http 


//www 


f eynarts 


de/f ormcalc 




http 


//www 


f eynarts 


de/looptools 
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